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Abstract. We propose a modified pairing functional for nuclear structure calculations which avoids the 
abrupt phase transition between pairing and non-pairing states. The intended application is the description 
of nuclear collective motion where the smoothing of the transition is compulsory to remove singularities. 
The stabilized pairing functional allows a thoroughly variational formulation, unlike the Lipkin-Nogami 
(LN) scheme which is often used for the purpose of smoothing. First applications to nuclear ground states 
and collective excitations prove the reliability and efficiency of the proposed stabilized pairing. 
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] 1 Introduction 

I Pairing is a key ingredient in a mean field description 
, of nuclear structure [1]. The actual handling has much 
' evolved. Early models referred to simple forms assuming a 
I constant gap or constant pairing matrix element [2] . Mod- 
■ ern nuclear mean field calculations use a fully fledged pair- 
ing force, often even density dependent, for an up to date 
. review see [3]. In any case, the theoretical description re- 
' lies on an extended notion of a mean field, now including 
, the two-quasi-particle density and gap-potential also at 
' mean field level. That allows to include a crucial part of 
' two-particle correlations at low expense. The price to be 
. paid for that simplification is the occurrence of a pairing 
' phase transition. Whenever the pairing strength crosses a 
, critical value, the system jumps from a pairing to a non- 
> pairing state or vice versa. For finite system such a sudden 
transition is an artefact [3]. A fully correlated treatment 
would produce correctly smoothed transitions. That goal 
of smoothing can already be achieved when working with 
a particle-number projected pairing state, provided the 
mean-field variation is done after projection (VAP). For 
then the system takes advantage of the pairing channel to 
incorporate some correlations "knowing" that the exact 
particle number will be restored. But the VAP method is 
cumbersome and not applicable in connection with en- 
ergy functionals [5]. Thus one employs usually simpler 
approaches which then are plagued by phase transitions. 
These are bearable when tracking ground states of nuclei 
where the change in nucleon number is discontinuous any- 
way. However, they turn into a severe problem when going 
along continuous changes as one does for large amplitude 
collective motion, see e.g. [3l[6j[7]. A sudden phase transi- 
tion along the collective evolution leads to unwanted sin- 
gularities in the description. Within the adiabatic treat- 
ment of large-amplitude collective motion, one encoun- 



ters, e.g., singularities in the collective masses. A phase 
transition would be even more disastrous for fully time- 
dependent mean-field calculations which enjoy a revival 
those days, see e.g. [8l|9]. Thus there is a need to smoothen 
the sudden pairing phase transition. A widely used recipe 
is the Lipkin-Nogami (LN) scheme pUUTT] which aims at 
an approximate treatment of particle-number projection 
and this approximation has the similar smoothing effects 
as the exact projection in VAP. Thus most of the mi- 
croscopic calculations of large-amplitude collective motion 
employ LN pairing to generate the underlying series of col- 
lectively deformed states, see e.g. [T2l[T3lfT4l[T5| . However, 
the LN prescription has still a few disadvantages in these 
large-scale applications: it introduces complicated terms 
and one more equation into the pairing scheme; the neces- 
sary feed-back on the self-consistent mean field [16] slows 
down the computational speed and degrades the numeri- 
cal stability; it is not strictly variational which makes the 
validity of the LN approach less controllable. Moreover, 
the feature of approximate particle-number projection is 
not really exploited. If a precise particle number is a mat- 
ter of concern, it is restored explicitely after the mean-field 
calculations. The superposition of mean-field states to a 
collective, or correlated, state requires anyway an explicit 
handling of particle number which is performed either by 
projection (e.g. |15| ) or by adjustment in the average (e.g. 
1141). As mentioned above, the formally perfect alterna- 
tive of particle-number projection with subsequent mean- 
field variation (VAP) is inhibited in connection with the 
Skyrme-Hartree-Fock functional There is thus a need 
for a conceptually simple, variationally consistent, and nu- 
merically stable recipe to avoid, or smoothen, the pairing 
phase transition, particularly for the simulation of large- 
amplitude collective motion. It is the aim of this paper 
to present and to discuss a stabilized pairing functional 
which avoids the phase transition while inducing a mini- 
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mum of changes to other nuclear properties. The idea is 
to modify directly the pairing-energy functional such that 
a vanishing pairing gap becomes associated with infinitely 
increasing energy while merging into the ordinary pairing- 
energy functional in regions of well developed pairing gap 
such that standard pairing properties are not affected. 
The stabilized scheme is illustrated on a schematic two- 
level model. First realistic application to nuclear struc- 
ture and collective excitations computed with the self- 
consistent Skyrmc-Hartree-Fock (SHF) method will also 
be presented. 



2 The stabilized functional 

The standard procedure starts from an energy functional, 
E ~ E^[ + -Epair, consisting out of a mean- field part 
and an pairing functional. The energy depends on a set 
of single-particle wavefunctions fk and occupation ampli- 
tudes Vk with complementing non-occupation amplitudes 
Uk = \/l — i^^. Variation with respect to (pi yields the 
mean-field equations and variation with respect to Vk the 
gap equation determining the pairing properties pilTT]. 
The gap equation does not always have finite solution. 
There can emerge a phase transition to zero pairing gap 
and zero pairing energy. This induces singularities in the 
derivatives which spoils many applications as, e.g., the de- 
scription of large-amplitude collective motion. The Lipkin- 
Nogami scheme (LN) avoids that sudden change by an ad- 
ditional term cx (iV^) in the energy, whereby N is the op- 
erator for particle number (for a discussion in connection 
with self-consistent mean- field models see [IS])- However, 
it has some disadvantages as discussed in the introduction. 
We propose to stabilize a finite pairing by introducing a 
counteracting term in the pairing energy functional. The 
stabilized pairing functional reads 
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It can be applied to any given pairing functional impair- 
It guarantees always a finite pairing because the term oc 
-^pair grows hugc whcu pairing is on the way to breakdown. 
The stabilized function llj introduces a new parameter, 
the cutoff pairing energy -Ecutp which has to be chosen to 
deliver sufficient stabilization with minimal effects in the 
well pairing regime. The stabilized equations are derived 
as before by variation of the total energy, now 
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This yields the stabilized coupled mean-field and gap equa- 
tions 
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where Ak is composed from the Uk,Vk as in regular BCS 
approach for the given functional. The detailed form de- 
pends on the actual system. It will be specified later on in 
connection with the examples. 



3 Test in a two- level model 

For a first exploration we employ a simple-most model for 
mean field plus pairing. We consider N particles in two 
energy levels, both iV-fold degenerated. The levels are la- 
beled by (s,to) where s E {u, 1} stands for the principle 
quantum number and m = ±1, ±N/2 accounts for the 
(degenerated) magnetic quantum numbers. The energy 
difference between upper and lower states, E^i = Eu ~ Si 
simulates the shell spacing. Pairing is added with a con- 
stant pairing force [5] . The occupation amplitudes are also 
degenerated Vsm = Vs- For symmetry reasons, we have 
vi = V, Vu = Vl — v'^, and u 
gethcr the total energy 



y/l — v'^. This makes alto- 
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The mean-field variation is obsolete because the single par- 
ticle energies eu,i are given. It remains the variation of 
pairing properties which yields the gap equation E^iuv + 



A{u^ -V' 
becomes 



where A = GNuv. Its solution for the gap 
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The decisive parameter for the model is the ratio NG/Eui- 
Pairing breaks down for NG/E^i < 1. The breakdown 
should now be hindered by introducing the stabilized func- 
tional ll]) which yields the stabilized gap as 
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This is a rather involved equation for the gap because -Epair 
on the right hand side depends implicitly on A. But the 
crucial result is that the expression under the square root 
can never vanish and thus the gap equation has always a 
finite real solution. 

Figure [T] shows results of the two-level model for phase 
space TV = 20. The pairing gaps are given in the upper 
panel. BCS produces a phase transition at relative cou- 
pling strength G/Eu\ w 0.05 and LN clearly removes that 
discontinuity producing a gap which smoothly approaches 
zero in the limit G — > 0. In the regime of active pair- 
ing, the LN gap approaches the BCS gap as it should. 
Both results from the stabilized functional ([1]) smoothen 
the transition, at first glance similar to LN. However, the 
functional behavior differs in detail. The trend at small G 
produces more efficient stabilization of the pairing gap, as 
can be seen, e.g., from the result for the extremely small 
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Fig. 1. Upper: Dependence of the pairing gap on the coupling 
strength for various methods in comparison, BCS, LN and the 
stabilized pairing with two different cutoff parameters as indi- 
cated. Lower: The total energy of the system as function of the 
coupling strength. 



G/Eyi = 0.01. At the side of active pairing, the stabi- 
lized functional also converges to the BCS curve where the 
speed of convergence improves with decreasing i?cutp- The 
goal is to have a slow decrease of A for G — > combined 
with a fast convergence to BCS for large G. The stabilized 
pairing promises to provide here a better compromise than 
LN. The lower panel of figure [T] shows the binding ener- 
gies. LN follows the BCS curve with a maximum deviation 
of 2% near the critical value G/E^i ~ 0.05. The stabi- 
lized pairing shows also deviation which, however, extend 
further below the critical value. The maximum deviation 
depends sensitively on the cutoff energy. The smaller cut- 
off produces an error similar to LN and the larger cutoff 
has a larger error. Thus a good compromise will dwell at 
the lower side of cutoff energies. The actual choice will 
depend, of course, on the actual system and the intended 
applications. 

The two-level model is based on an explicit Hamilto- 
nian. This allows to perform particle-number projection. 
The most advanced approach is mean-field variation af- 
ter projection (VAP). We compare that with the results 
from mean-field calculations where projection is done af- 
ter mean-field variation. Figure [2] shows the results drawn 
as relative deviation from VAP because the differences 
for the projected results are too small to be seen when 
showing the total energies as such. The BCS case devi- 
ates most strongly below the critical coupling strength G. 
Pairing was not active here and no energy can be gained 
from pairing with subsequent projection. The deviation 
stays at the level of 3-4 % for pure mean- field results. 
The unprojected BCS curve shows the typical size and 
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Fig. 2. Comparison of various approaches with variation af- 
ter particle-number projection (VAP) for the total energy as 
observable drawn as % deviation from the VAP result ver- 
sus the normalized coupling strength. Test case is again a 
system with particle number = 20. The four cases, BCS, 
LN and stabilized BCS at two different cutoff energies employ 
particle-number projected calculations applied to the given 
mean field state (projection after variation = PAV). An un- 
projected mean- field result is also given for the case of BCS. 



trend which looks very similar also for the (unprojected) 
other methods. The projected BCS result converges above 
the critical strength rapidly to a very small deviation. LN 
and stabilized gap show very small deviations throughout. 
The figure confirms that LN with subsequcnz projection 
is an efficient approach to full VAP. And it shows that 
the newly developed stabilized gap performs equally well, 
and sometimes better, in that perspective. The case of 
EcMtp/Eiii = 0.6 shows that proper tuning may even en- 
hance the quality. On the other hand, being satisfied with 
LN quality means that one disposes of a broad band of 
choices for -Ecutp/i'ui which one can use to optimize other 
aspects as, e.g., stability. 



4 Nuclear ground states and collective 
dynamics 

4.1 The modified functional and BCS equations 

Self-consistent nuclear mean-field models are based on ef- 
fective energy-density functionals. The three most widely 
used models are the relativistic mean-field model (for re- 
views see e.g. p^[T9l[20] ). the Gogny force (see e.g. [2T]). 
and the Skyrme-Hartree-Fock (SHF) approach (for a re- 
cent review see e.g. [3]). We concentrate here on SHF. 
Although there had been attempts to incorporate pairing 
directly into the Skyrme force [22], most functionals keep 
mean-field and pairing properties separate. Recent SHF 
calculations employ generally a pairing functional derived 
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from a zero-range pairing force p3l[24] . It reads 



> 

0) 



0.5 



^P- = iE / d'rV^i^)xl{r) , (6a) 

(6b) 



Xq(i") = 2^/feUfcUfc|<^fe(r 

fe>0 



where fk is a smooth eutofF function which hmits the pair- 
ing space to an energy region up to 5 MeV above the Fermi 
energy [31125] . That pairing functional is added to the SHF 
mean-field energy functional which is rather lengthy. We 
refer to |3l[26j for details. The state-dependent pairing gap 
from this functional becomes 
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The gap equations derived from that functional show also 
the possible breakdown of pairing. We modify the func- 
tional ([6a| to the stabilized pairing ([1]). As a consequence, 
we deal now with the stabilized equations ([2]). These will 
be solved in spherical symmetry for the ground states of 
spherical nuclei [27j and in axial symmetry for computing 
the collective deformation path and subsequent collective 
dynamics [M]. There is a broad variety of parameteri- 
zations of the SHF mean-field functional [2- The actual 
choice is not crucial for the present tests of the stabilized 
pairing functional. We use here the parameterization SkI3 
[28]. 

The pairing gap from the functional ([Sa]) becomes state 
dependent as A^- In order to have one number character- 
izing the gap, we define an average pairing gap as 



J2keq fkUkVkAk 

Y.keq fkUkVk 



(7) 



The weight factor UkVk concentrates averaging to a region 
near the Fermi energy where pairing is most relevant. That 
definition is found to be the most appropriate for compar- 
isons 1^. 



4.2 Results and discussion 

The uppermost two panels of figure [3] show proton and 
neutron pairing gaps along the semi-magic chain of Sn 
isotopes. The proton number Z=50 is magic throughout 
the isotopic chain. Thus BCS proton pairing breaks down 
for all isotopes while LN and stabilized pairing produce 
a small but finite gap. The actual size of the stabilized 
gap depends somewhat on the chosen cutoff whereby the 
lower i?cutp produces results very close to LN. BCS neu- 
tron pairing drops to zero for the two magic neutron num- 
bers iV = 50 and = 82, has substantial pairing mid 
shell, and a short transitional region near the magic num- 
bers. LN produces, of course, a finite gap at the magic 
points and tries to come close to BCS when going away 
from the shell closures. The results from stabilized pair- 
ing do also produce finite gaps at the shell closures. The 
actual value of the gaps there depend again somewhat on 
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Fig. 3. Ground state properties for the chain of Sn isotopes 
computed with BCS, LN and stabilized pairing at two -Ecutp. 
Upper two panels: Average proton and neutron pairing gaps. 
Middle panel: difference 5Eb of binding energy from stabilized 
pairing or LN to the BCS binding energy. Lowest panel: differ- 
ence Srin-is for the charge r.m.s. radius. 



^'cutp- Outside the magic points, the gaps from stabilized 
pairing come very close to the BCS values, for the small 
as well as for the larger cutoff, and both closer than LN. 

The lower two panels in figure [3] show the effect on 
bulk observables. Results are drawn as difference to the 
BCS results in order to fit the small overall changes into a 
graphical representation. The changes on the total bind- 
ing energy (middle panel) have an offset along all isotopes 
and a peak at the magic neutron numbers. Offset and 
peak are expected because the proton stabilization is al- 
ways present and neutron stabilization adds to the effect 
for N=50 and N=82. The average size of the modification 
is close to the LN values for the smaller cutoff. We have to 
mention, however, that the energies from LN are an uncer- 
tain quantity because the LN scheme is not variationally 
consistent. It is thus not fully clear what really should be 
plotted here. The stabilized pairing, on the other hand, 
has by construction a well defined energy. And a change 
of the order of 0.2-0.5 MeV is well acceptable in view of a 
typical uncertainty of 0.7-1 McV for Skyrmc forces. The 
effect on radii (lowest panel) is extremely small as com- 
pared to typical uncertainties of about 0.02-0.04 fm [3]. 
The same holds for the whole density and deduced ob- 
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Fig. 4. Proton and neutron gaps (upper two), collective mass 
for motion along 020 (middle), and difference of the collec- 
tive potential from stabilized pairing or LN to that from BCS 
(lower) along the quadrupole deformation path for ^^^Sn com- 
puted with BCS, LN and stabilized pairing. 



servables as difFraction radius and surface thickness 
[30]. 

Low lying nuclear 2+ states are related to large am- 
plitude collective motion. The nucleus vibrates along a 
collective deformation path which consists out of a con- 
tinuous series of quadrupole deformed states. The path 
is generated by adding a quadrupole constraint to the 
mean field equations sampling the various quadrupole de- 
formations; from the given path, one can derive then the 
ingredients for the collective dynamics, the potential as 
expectation value of the energy and the collective mass 
from double commutators with the dynamical boost, for 
details see, e.g. [7lfT4]. Figure |4] shows results along the 
axially symmetric quadrupole deformation path of ^^^Sn. 
BCS shows the breakdown of neutron pairing near the 
spherical shapes and the mass jumps immediately to large 
values because a pure Slater determinant is more resistant 
to quadrupole motion than a well pairing BCS state. The 
sudden changes make the BCS mass unusable for further 
processing in the collective Hamiltonian. LN and the sta- 
bilized functionals succeed in preventing the breakdown 
of pairing and in producing a smooth transition. They be- 
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Fig. 5. Systematic of quadrupole-coUective properties along 
the chain of Sn isotopes computed with LN and stabilized pair- 
ing at two -Ecutp. Upper: The low- lying i5(2i-|-) excitation en- 
ergies. Lower: Collective correlation energies. 



have similar what the collective mass and the changes in 
energy is concerned. It is to be noted that the stabilized 
pairing remains closer to the BCS mass in the regime of 
shell closure. The net effect of all these small differences 
on the collective excitations will be an average over all 
deformations. 

Figure [4] has shown the impact of the stabilized pair- 
ing functional and its -Ecutp on the constituents of the col- 
lective Hamiltonian, mass and potential. Figure [5] shows 
the net effect on the final states. The results for LN and 
stabilized pairing are again similar in general with a few 
interesting differences in detail. The sudden reduction of 
the E{2'^) when departing from the magic iV = 82 is very 
badly reproduced by LN which produces a much smoothed 
transition. The stabilized pairing performs better in that 
respect. The mid-shell level is reached already at = 78 
and = 80 also does a larger down-step. The trend is 
the better reproduced the lower -Ecutp- The excitation en- 
ergies at the magic point and mid-shell are similar for all 
methods which proves that the modified pairing recipes 
do not affect too much the physics of the collective exci- 
tations. The sometimes larger effects seen in the details 
of collective mass and potential (figure 2]) seem to average 
out in the final collective energies. The results for the cor- 
relation energies are even more robust. The values differ 
in tiny detail, but that is not significant because the cor- 
relation energy is anyway a small correction to the total 
energy. 
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The most important result coneerns a praetieal aspeet, 
not shown in the above figures: Stabilized pairing is sim- 
pler to implement numerically and it produces a faster 
and much more stable iteration. This amounts to a fac- 
tor three faster computation of the collective deformation 
paths which is a substantial gain for large scale applica- 
tions. 



5 Conclusions 



We have investigated a new scheme to smoothen the pair- 
ing phase transition in nuclei. The stabilized pairing is 
achieved by modifying the pairing functional such that 
the pairing energy diverges for vanishing gap. The modifi- 
cation is restricted to the regime of small coupling and has 
little effect elsewise. The intended effect should come out 
similar to what is usually achieved by the Lipkin-Nogami 
(LN) scheme. 

The operation of the stabilized functional was demon- 
strated in a schematic two-level model and compared to 
the LN. It was shown that a smooth transition to vanish- 
ing gap can indeed be achieved by the stabilized pairing 
functional while the convergence towards BCS results in 
the well pairing regime is even somewhat faster for stabi- 
lized pairing as compared to LN. The BCS ground state 
with stabilized pairing can also serve very well as starting 
point for subsequent particle number projection. The re- 
sults come very close to fully fledged mean-field variation 
after projection. 

The stabilized functional has been applied to a bunch 
of realistic nuclear structure calculations, ground states of 
spherical nuclei and low-energy collective dynamics. All 
results show that stabilized pairing is a valid alternative to 
LN. The effects are very similar in general. There are some 
small differences in detail which give a slight preference to 
stabilized pairing. 

The results depend somewhat on the cutoff parameter 
inherent in the model. Pairing properties are, of course, 
more sensitive, but bulk properties, fortunately, not so 
much. A good choice which has small side-effecs and comes 
in many aspects close to LN is a small value around i?cutp = 
O.SMeV. 

The major advantages of stabilized pairing lie in a for- 
mal and in a numerical feature. The formal advantage is 
a thoroughly variational formulation. The modification is 
done only at the side of the functional which provides a 
clean definition of the energy. The numerical gain is dra- 
matic. Stabilized pairing provides factor three faster con- 
vergence than with LN due to the simpler algorithm pos- 
sible (remaining practically identical to the robust BCS 
scheme). This is the strongest argument in favor of the 
new scheme. A stabilized pairing functional with fully vari- 
ational formulation will become inevitable in future time- 
dependent applications where a possible phase transition 
would be disastrous and where LN is not applicable. 
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